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ABSTRACT 

Stellar population models of absorption line indices are an important tool for the 
analysis of stellar population spectra. They are most accurately modelled through 
empirical calibrations of absorption line indices with the stellar parameters effective 
temperature, metallicity, and surface gravity, the so-called fitting functions. Here we 
present new empirical fitting functions for the 25 optical Lick absorption line indices 
based on the new stellar library MILES. The major improvements with respect to the 
Lick/IDS library are the better sampling of stellar parameter space, a generally higher 
signal- to-noise, and a careful flux calibration. In fact we find that errors on individual 
index measurements in MILES are considerably smaller than in Lick/IDS. Instead we 
find the rms of the residuals between the final fitting functions and the data to be 
dominated by errors in the stellar parameters. We provide fitting functions for both 
Lick/IDS and MILES spectral resolutions, and compare our results with other fitting 
functions in the literature. A Fortran 90 code is available online in order to simplify 
the implementation in stellar population models. We further calculate the offsets in 
index measurements between the Lick/IDS system to a flux calibrated system. For 
this purpose we use the three libraries MILES, ELODIE, and STELIB. We flnd that 
offsets are negligible in some cases, most notably for the widely used indices H/3, 
Mgb, Fe5270, and Fe5335. In a number of cases, however, the difference between flux 
calibrated library and Lick/IDS is signiflcant with the offsets depending on index 
strengths. Interestingly, there is no general agreement between the three libraries for 
a large number of indices, which hampers the derivation of a universal offset between 
the Lick/IDS and flux calibrated systems. 
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1 INTRODUCTION 

Stellar population models of absorption line indices 
are a key tool for the analysis of star cluster and 
galaxy absorption spectra. They are used to derive 
the fundamental stellar population properties such as 
age, metallicity and element abundance ratios. In par- 
ticular, optical absorption line diagnostics in the spec- 
tra of evolved stellar populations have successfully been 
adopted in the past in studies on g alaxy evolution (e.g. 
Worthev. Faber fc Gonz alez 1992; Davies, Sa dler fc Peletieil 
19971 



19931: IVazdekis et al 



Worthevlll998l: iTrager et al 



jKuntschncr^ 
199 8: Hcnrv ' 



DaviesI 



Worthe 



1998 



199S 



KuntschneJ I2OO0I : iTrager et al.l I200Q : iThomas et al. 2005 
and globular clus t er formation (e.g. iKissler-Patig 1998 



2005 



Forbes et all I2OOII : iKuntschnej I2OO2I : iBrodie et al] 
Puzia et al.l boOSi). The Lick/IDS system (iBurstein et all 
19841 : iFaber et al.lll985l: iGorgas et al.lll993l: (worthev et all 
1994 : Iworthev fc Ottavianilll997l : ITrager et al.lll998h is the 



standard set of absorption line indices that has been used 
extensively during the last two decades for studying absorp- 
tion features of stellar populations. This system consists of 
index definitions for 25 prominent absorption features be- 
tween 4000 and 6500 A present in the spectra of evolved 
stellar populations. 

For studies of galaxy and star cluster evolution, 
absorption lines need to be modelled for stellar popula- 



tions (e.g. 


MarastonI Il998l. l2005l: 


Bruzual & CharlotI l2003l: 


Vazdekis 


I999I: Leitherer et al 




19991: Worthev et all 


19941: 1 Worthev & Ottavianil Il997 


: iRose et al.l 


199J; 
20031; 


Trager et all I2OO0I: iThomas. Maraston & Bended 


Thomas. Maraston fc KornI 2004) 


A convenient way goes 



through the use of empirical calibrations. This is motivated 
by the fact that theoretical model atmospheres are known to 
suffer fro m incomplete line lists and continuum uncertain- 



ties, (e.g. Korn. Maraston fc Thomas! 20051: Coelho et al.l 



20071 : iRodn'guez- Merino et al . 2005; Le e. Worthev fc Potter! 
2OO9I : IWalcher et al.l bOol T Empirical calibrations on the 
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other hand have the disadvantage to be hardwired to the 
chemical abundance pattern of the Milky Way, which can 
be overcome in a semi-empirical approac h as in the rn odels 
by iTrager et all (|2000l) . I Thomas et al] (|2003l . |2004| ) and 
ISchiavonI {20m 

Empirical calibrations can be inserted in the models 
in two ways. In the first and most widely used approach, 
absorption line indices enter stellar population modelling 
through calibrations of the empirical relationship between 
the indices and the stellar atmospheric parameters Tcff, 
logg and [Fe/H] as provided by stellar libraries. As these 
calibrations are usually obtained through polynomial 
fitting procedures they are commonly referred to as 'fitting 
functions'. The quality of the final stellar population 
model critically depends on the accuracy with which these 
relationships can be inferred from stellar libraries, i.e. the 
coverage of stellar parameter space and the reliability of 
the index measurements. The computational procedure 
with which the fitting functions are determined is a further 
crucial step in producing accurate models. A number of 
studies in the literature are devoted to such empirical 
calibrations for various stellar libraries, either for the Lick 
indices, parts of t he Lick indices or other prominent ab- 
sorption features jBuzzoni. Garib oldi fc Mantegazza 19921: 



Buzzoni. Mantegazza fc Garibold i 1994i; IWorthev et al 



1994' 



Borgcs ct al. 1995 



i'; 'Gorgas ct al. 19991: ICenarro et al 



2002 : Schiavo n 2007: Maraston et al. 2008 

Alternatively to the use of fitting functions, absorp- 
tion line indices can be measured directly on the synthetic 
spectral energy distribution (SED) from stellar population 
models that are based on empirical stellar libraries. The 
benefit of this method is that the full SED can be com- 
pared pixel-by-pixe l to observations (e.g. IPanter et al ] |2007l : 
iToieiro et alfi2007l ). 

The major strength of fitting functions, instead, lies in 
the fact that they allow for interpolation between well pop- 
ulated regions of stellar parameter space which increases 
the accuracy of the model in stellar parameter space that 
is only sparsely sampled by empirical stellar libraries. More- 
over, each absorption index or spectral feature is represented 
by an individual fitting function, which is optimised to best 
reproduce its behaviour in stellar parameter space. Fitting 
functions are also easier to implement in a stellar popula- 
tion synthesis code, and models based on fitting functions 

are better comparable. 

The w idely used fitting functions of IWorthev et al.l 

l|l994h and IWorthev fc Ottavianil (Il997l) are based on the 



Lick/ IDS stellar library ( Burstein et al.| [l984: Fa ber et al.l 
119851). They are adopted in most stellar population mod- 
els llWorthevlll994l: IVa^dekis et al.lll996l: iTrager et al.ll2000l : 



iThomas et al.ll2003l . l2004l2005l : lAnnibah et al.ll200'i1) in the 

literature. Other fitting func tions based on th e same stellar 
hbrary exist (B uzzoni et al.1 ,1992. 1994; Borg cs et aLlliggSl ) 
and lead to overal l consistent results in the final stellar 
population model (|Maraston et al.l |2003| ). Major progress 
has been made with the advent of a new generat i on of 
stellar libraries djoncs '1999*; 'Prugnicl fc Soubiranl [iooH : 



iLe Borgne et al.] [2003; Sanchcz-Blazqucz ct al. 2006|) that 
have led to considerable improvements regarding coverage of 
stellar parameter space, spectral resolution, signal-to-noise 
ratio, and flux calibration. 

In particular the latter is a critical step forward. As the 



Lick/IDS system is not flux calibrated, observations have 
to be re-calibrated onto the Lick/IDS system through com- 
parison with Lick standard stars. This requirement hampers 
the analysis of data samples for which spectra of such cali- 
bration stars are either not available at sufficient quality or 
do not cover the appropriate rest-frame wavelength range. 
This problem is most imminent in high redshift observations 
and in g alaxy redshift su rveys such as the Sloan Digital Sky 
burvey (|York et al.ll2000l ). The new flux calibrated libraries 
allow the analysis of flux calibrated spectra at any redshift 
without spectroscopic standard stars. 

Flux-calibrated stellar libraries in the literature that 
are suitabl e for stellar population modelling include 
the Jones l|Jonesl [1999) ■ EL OP IE (IPrugniel fc Soubiranl 
l200j). STELIB (iLe Borgne et af] |2003| ) and MILES 
ijSanchez-Blazguez et al.ll2006l ) libraries. The MILES library 
is particularly well suited for stellar population modelling 
of absorption line indices owing to its favourable combina- 
tion of spectral resolution, wavelength range, stellar param- 
eter coverage, and quality of fiux calibration. In this paper 
we present new Lick index fitting functions based on the 
MILES stellar library. To take advantage of the full spectral 
resolution of the MILES library we have produced fitting 
functions for both the lower Lick/IDS resolution (8 — 11 A 
FWHM) and the higher resolution of the MILES library (2.3 
A FWHM). A new version of the TMB stellar population 
model of absorption line indices based on these new fitting 
functions will be presented in a subsequent paper. 

This paper is organized as follows. In Section [2] we 
present the Lick indices measured on the MILES library and 
a quality evaluation of the index measurements. We discuss 
offsets between the flux calibrated MILES and the Lick/IDS 
systems. The empirical fitting method is presented in Sec- 
tion|3]along with the resulting fitting functions. In Section|4] 
we compare the fitting functions of this work with fitting 
functions from the literature. We summarise in Section [S] 



2 THE MILES STELLAR LIBRARY 



The MILES hbrary ([Sanchez-Blazguez et al.ll2006l ) consists 
of 985 stars with spectra in a wavelength range of 3525- 
7500 A, well covering the Li ck indices, and with a spec- 
tral resolution of 2.3 A (see ISanchez-Blazguez et al.l |2006| 
for further details). Important for the aim of this work 
is the careful flux-calibration of the MILES spectra. Also, 
ISanchez-Blazguez et al.l (|2006l ) selected the sample of stars 
to fill the gaps in stellar parameter space covered by previous 
stellar libraries. This makes the MILES library particularly 
suitable for modelling absorption line indices of stellar pop- 
ulations. 

Stellar parameter estimates in the literature show a 
scatter due to varyin g methods applied, as discussed in 
iMaraston et al. [ l|2003l ) for [Fe/H]. The steUar parameters 
(Tefi, logg o r [Fe/H]) for the star s in the MILES library are 
presented in lCenarro et al.l(|2007D . where estimates from the 
literature have been used and put on a homogeneous scale. 3 
of 985 stars have no available estimates for none of the stel- 
lar parameters Teff, logg or [Fe/H]. 35 stars lack estimates 
only for [Fe/H] and are located in sparsely populated regions 
at the ends of the Tcff range. The stars have therefore been 
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Table 1. Typical Lick index errors and offsets to the Lick/IDS library. M-it and L-u corresponds to index errors at the resolution of the 
MILES and Lick/IDS libraries, respectively. T98-(t are the index errors presented in lTrager et al.l lll99Sl ) for the Lick/IDS library. Im, 
are indices measured on the libraries {MILES, ELODIE and STELIB) for which offsets to the Lick/IDS library are presented. Il are 
indices measured on the Lick/IDS library. 



INDEX Error Offset 1^, = a ■ 1^^^^ + b 

MILES ELODIE STELIB 



i 


NAME 


M-cr 




L-a 


T98-0- 


a 


b 


a 


b 




a 


b 


1 


KSa 


0.164 





125 


0.64 


0.960 


-0.054 


0.955 


0.721 





940 


0.823 


2 


HSf 


0.093 





075 


0.40 


0.965 


0.049 


0.936 


0.397 





956 


0.242 


3 


CNi 


0.0042 





0038 


0.018 


0.912 


0.008 


0.897 


-0.012 





986 


-0.010 


4 


CN2 


0.0050 





0042 


0.019 


0.907 


0.006 


0.900 


-0.008 





985 


-0.013 


5 


Ca4227 


0.063 





047 


0.25 


0.904 


0.074 


0.771 


0.163 





918 


-0.057 


6 


G4300 


0.112 





093 


0.33 


0.858 


0.625 


0.870 


0.646 





924 


0.565 


7 


H7A 


0.142 





107 


0.48 


0.976 


-0.148 


0.967 


-0.057 


1 


022 


-0.735 


8 


H7f 


0.069 





059 


0.33 


0.963 


-0.038 


0.962 


0.016 





999 


-0.238 


9 


Fe4383 


0.155 





127 


0.46 


0.932 


-0.220 


0.929 


-0.184 





915 


0.796 


10 


Ca4455 


0.073 





056 


0.22 


0.747 


-0.067 


0.785 


-0.105 





891 


-0.228 


11 


Fe4531 


0.122 





096 


0.37 


0.857 


0.290 


0.838 


0.390 





877 


-0.002 


12 


C24668 


0.179 





156 


0.57 


0.903 


0.484 


0.913 


0.295 





992 


0.512 


13 


H/3 


0.063 





051 


0.19 


0.981 


0.126 


0.996 


0.015 


1 


004 


0.032 


14 


Fe5015 


0.139 





115 


0.41 


0.902 


0.084 


0.926 


0.178 





989 


0.168 


15 


Mgi 


0.0017 





0013 


0.006 


0.911 


0.0004 


0.923 


0.005 





903 


-0.009 


16 


Mg2 


0.0023 





0014 


0.007 


0.918 


-0.003 


0.940 


0.0006 





960 


-0.013 


17 


Mgb 


0.053 





045 


0.20 


0.964 


0.108 


0.935 


0.247 


1 


003 


-0.026 


18 


Fe5270 


0.058 





047 


0.24 


0.923 


0.101 


0.919 


0.180 





932 


0.173 


19 


Fe5335 


0.063 





044 


0.22 


0.960 


0.135 


0.963 


0.032 





946 


0.110 


20 


Fe5406 


0.044 





031 


0.18 


0.874 


0.269 


0.913 


0.165 





853 


0.264 


21 


Fe5709 


0.060 





050 


0.16 


0.979 


-0.026 


0.907 


0.015 


1 


019 


-0.046 


22 


Fe5782 


0.057 





043 


0.19 


0.920 


0.037 


0.879 


-0.004 





906 


0.088 


23 


Na D 


0.082 





064 


0.21 


0.990 


-0.162 


0.979 


-0.069 





993 


-0.071 


24 


TiOi 


0.0021 





0017 


0.006 


0.918 


-0.005 


0.895 


-0.006 





918 


0.0003 


25 


TiOa 


0.0022 





0016 


0.006 


0.904 


0.0007 


0.912 


0.005 





940 


0.009 



assigned a solar metallicity to increase the number of data 
points. 



2.1 Empirical stellar Lick indices 

Our aim was to produce fitting functions both for the reso- 
lution of the MILES library (2.3 A) and for the resolution of 
the Lick/IDS library (8-11 A). We have therefore measured 
the 25 Lick indices directly on the original stellar spectra 
and on the spectra downgraded to t he Lick/IDS resolution 
described by the curve presented in IWorthev &: Ottavianil 



1997 ) . We have used t he index definitions fromlTrager et al.l 
199^ ) and also from IWorthev fc Ottavianil 

mm for the 



higher order Balmer lines (RSa, HSf, K'Ja and H7f). Ob- 
servational errors and offsets to the Lick/IDS library are 
described in the following paragraphs. 



2.1.1 Observational index errors 

We have derived typical observational index errors in order 
to evaluate the quality of our index measurements. To this 
end we have used pixel l-a observational errors (P. Sanchez- 
Blaszquez private communication) to perturb each stellar 
spectrum, both at MILES resolution and Lick/IDS resolu- 
tion, through 600 Monte Carlo realizations. We have then 
measured the 25 Lick indices for each perturbed spectrum 
and determined l-a errors for each index by using the spread 



in index measurements from the realizations. The index er- 
rors of the individual stellar spectra are used for weighting 
the least square fits when deriving both the ofi^sets to the 
Lick system (Section I2.1.2p and the fitting functions (Sec- 
tion [3)1. 

Trends between the index errors and the atmospheric 
parameters or line-strength indices can in principle bias the 
fits, but we have found such trends not to affect the re- 
sults. Only for the Balmer indices we find weak trends of 
increasing errors with decreasing temperature and decreas- 
ing index strength. No trends with logg and Fe/H are found 
for the Balmer indices. These weak trends can probably be 
explained with higher S/N for bright hot stars where the 
Balmer indices increase significantly in strength. Since we 
compute the fitting functions in bins of temperature, these 
trends have no significant effects on the final fitting func- 
tions. 

The final l-a typical index errors were determined by 
taking the median error of the whole stellar library for each 
index. The typical index errors are presented in Table[T]both 
for MILES and Lick/IDS resolution. Compare d to the typi- 
cal in dex errors for the Lick/IDS stellar library l|Trager et al.l 
'19981), also included in Table [TJ we find the errors of the 
MILES library to have improved significantly. The stars of 
the Lick/IDS library were observed about thirty years before 
the MILES library. Considering the technical development 
in thirty years time, an improvement in the measured indices 
ought to be expected. 





Figure 1. Index by index comparison between index strengths measured on the MILES hbrary and the Lick/IDS library. Each 
panel shows the residual as a function of index strength. Included are also least-square fits of the residuals (black lines, coefficients 
presented in Table [TJ that show clear index strength dependent offsets between the two libraries. Red crosses are sigma clipped 
data points in the least-square fitting routine. Typical index errors are indicated in the right top corners (see text for more details). 
Included are also offsets to the Lick/IDS library derived in this work fo r the ELODIE (magenta lines) and STELIB (green lines) 
libraries, as well as the offsets derived in lSanchez-Blazguez et al.l ||2009| ) for the MILES library (dotted lines). 



2.1.2 Lick Index offsets 

We have computed Lick index offsets between the MILES 
library and the Lick/IDS library using the stars in com- 
mon between the two libraries. These offsets can be 
used for comparisons between models based o n this work 
with models based on th e IWorthev et al.l () 19941 ) and 
IWorthev fc Ottavianil (|l997l ) fitting functions. The offsets 
are also used in Section |4]to compare the fitting funct i ons o f 
this work with the fitting functio ns of IWorthev et all (|l994l ) 
and lWorthev fc Qttaviajiil ([l993). 

Fig. [T] shows index by index comparisons for the resid- 
uals between the index meas urements of the two l i brarie s 
as function o f inde x stre ngth. Worthey fc Ottavianil l|l997h , 



iKuntschneil (|200ll ) and ISchiavonI (f200'iD computed zero- 
point offsets to the Lick/IDS library, while IPuzia et al.l 
l|2002 ) computed their offsets as 2nd order least-square fits. 
For most indices we find index strength dependent residuals 
between the two libraries (Fig. We have therefore com- 
puted the offsets using a sigma-clipping linear least-square 
fitting routine, weighted with the individual index errors de- 
rived in Section [2.1.11 The slope and intercept of these fits 



are presented in Table [T] and also included in Fig. [T] (black 
solid lines). Sigma-clipped data points are indicated with red 
crosses in Fig. [T] and the error bars are the 1-a index errors 
presented in Section 12.1.11 The error bars along the x-axes 
are represented by the index errors derived for the MILES 
library, while the error bars along the y-axes are represented 
by the combined errors of the MILES and Lick/IDS libraries 
in quadrature. 

Extreme outliers, i.e. data points that clearly showed 
strong deviating values compared to the bulk of data points, 
were removed prior to running the fitting routine, in order 
to avoid stars with anomalous index strengths to affect the 
final fits. For three indices (C24668, Fe5015 and Mgb) we 
found offsets at particularly high index strengths that devi- 
ated from the offset trends for the majority of data points. 
The low number of data points at these index strengths and 
the absence of data points at intermediate index strengths 
induced a bias in the derived offsets. The data points at par- 
ticularly high index strengths were therefore discarded when 
determining the final offsets. 

Offset s between the MILES and the Lick/IDS library 
derived in ISanchez-Blazguez et ahl (|2009l ) are also included 
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Figure 1 - continued 



in Fig. [T] (dotted lines). These offsets and the offsets derived 
in this work for the MILES library show in general very good 
agreement. Differences greater than the l-o index errors are 
mainly found for Mgb, one of the indices for which we ex- 
cluded data points at high index strengths due to deviations 
in offset trends. Noticeable offset differences, but still within 
the l-cr index errors, are also found for Fe5015, Fe5709, NaD, 
TiOi and Ti02. Only small deviations are found between the 
offsets derived in lSanchez-Blazauez et al] (|2009l l and in this 
work for the rest of the indices. 

For comparison, and for the derivation of a possi- 
ble universal offset between flux-calibrated system and the 
Lick/IDS system, we have also determined offsets to the 
Lick/IDS library for two other flux-calibrated stell ar li- 
braries, n amely ELODIE jPrugniel fc Soubiranll2001^ and 
STELIB (|Le Borgne et al.n2003l ). These offsets were deter- 



mined using the same procedure as described above for the 
MILES library, except that no individual index errors were 
used as weights in the least-square fitting. For STELIB the 
lack of information did not allow for a computation of in- 
dex errors, while the derived index errors for ELODIE were 
found to be unreliable as they showed unrealisticly small 
values. Since we only found small deviations in the offsets 
derived for the MILES library when not weighting as com- 
pared to weighting the least-square fits, we compare the off- 
sets derived for all three libraries. 

In accordance with the MILES library we found in- 
dex strength dependent offsets also for the ELODIE and 
STELIB libraries. We found deviating offsets trends at high 
index values for the same indices as for the MILES library 
(C24668, Fe5015 and Mgb). 

The offsets derived for the ELODIE and STELIB li- 
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braries are also presented in Table[T]and Fig. [T] (magenta and 
green lines, respectively). Clearly, deviations in the offsets 
are found between the libraries, especially for the STELIB 
library compared to the other two libraries. However, the 
STELIB library is also the library having the least number 
of stars in common with the Lick/IDS library, giving a higher 
statistical uncertainty in the derived offsets. The STELIB 
library only has 44 stars in common with the Lick/IDS li- 
brary, while EL O DIE has 112 stars and the MILES library 
has 237 stars in common with the Lick/IDS library. Ca4227 
showed particularly strange behaviour with index strength 
and the accuracy of the final offsets for this index could be 
questionable. 

In Fig. [T] we find agreements within the l-cr index er- 
rors between the offsets derived for all three libraries for 
H/3, Mgb, Fe5270, Fe5335, Fe5406, Fe5709, Fe5782 and NaD. 
This implies a better agreement between all libraries at 
wavelengths redder than ~4800 A, with the exception for the 
broader molecular indices Mgi, Mg2, TiOi and Ti02 that 
show differences greater than the 1-a index errors, which is 
also found for Fe5015. Agreements between offsets derived 
for MILES and ELODIE only, well within the l-cr index 
errors, are found for G4300, H7A, H7F, Fe4383, Ca4455, 
Fe4531 and C24668. This instead implies a worse agreement 
between MILES and ELODIE at wavelengths bluer than 
~4250 A (H(5a, H(5f, CNi, CN2 and Ca4227), where we in 
general find inconsistencies between all three libraries. The 
significant deviation in offset between the libraries for sev- 
eral indices hamper the derivation of a universal offset be- 
tween the Lick/IDS and flux-calibrated systems as described 
by these libraries. 

This conclusion get s furth er support from the study of 
ISanchez-Blazguez et al.l (|2006l ) who show that offsets exist 
between the three flux-calibrated libraries MILES, STELIB 
and ELODIE. These offsets are generally in good agreement 
with the individual Lick offsets found in this work. 



3 FITTING FUNCTIONS 

In order to produce empirical fitting functions for the MILES 
library, we combine our measured Lick indices with the cor- 
responding stellar atmospheric parameters (see Section 
It is a complex task to find the best relationship between in- 
dices and stellar atmospheric parameters, with several meth- 
ods available in the literature. The method adopted in this 
work is presented in this section along with the derived fit- 
ting functions. 

A user friendly Fortran 90 code is available online at 
www.icg.port.ac.uk/~johanssj to make the implementation 
of our fitting functions easier. 



3.1 Fitting method 

The relationship between Lick index strengths and stellar 
parameters shows a complex behaviour, making it difficult 
to find one reliable empirical fitting function for the whole 
parameter space. To solve this problem the parameter space 
must be divided into subregions where local fitting func- 
tions can be computed. However, it is desirable to find the 
simplest set of fitting functions and achieve a final represen- 
tation of the data that is as accurate as possible. Hence 



the limits of the subregions have to be carefully chosen. 
It has also to be assured that adjacent subregions overlap, 
making smooth transitions possible. For these transitions 
we have adopted co sine-weighted interpolations following 
ICenarro et ail l|2002l ). The choice of subregions are discussed 
in Section [3.21 

Following the extensive numb er of publish e d fit- 



ting functions in the literature (iWorthev et ak. 
Gorgas et al ] I1999I: ICenarro et all I2OO2I : ISchiaronl 



1994; 



2007 



Maraston et al.l 20091 ) . we use a linear least square fitting 
routine to determine the local relationships as polynomials 
in the following way 



I {6, [Fe/H], logg) = ^ ' ' [Fe/H]'' ■ logg' 



(1) 



where j, k,l ^ and the atmospheric effective tem- 
perature is represented hy — 5040/Tcff. The represen- 
tation of Tcft using is chosen due to the wide range of 
spectral types in the stellar library. The number of terms 
in Eq. [T] can be made arbitrarily high. However, the goal 
is to find the best compromise between simplicity and ac- 
curacy by discarding terms with higher order polynomials 
that are negligible or induce unphysical behaviours. To this 
end several me t hods have been developed in the literature. 
IWorthev et al.l (|l994f ) presented a method to find the con- 
verging rms scatter by successively including terms and test 
if th e rms scatter w as significa ntly reduced by mean s of a F- 
test. lGorgas et al.1 (jl999i ') and ICenarro et all (|2002l ') instead 
test if each term signifi cantly differed from zero through a 
T-test. ISchiavonI (|2007l ') point out that both methods men- 
tioned above are sensitive to th e coverage of parameter 
space. Therefore ISchiavonI (|2007l ) combine the two meth- 
ods by first successively removing statistically insignificant 
terms and then interactively testing the remaining terms for 
unphysical behaviours and their effect on the rms scatter. 

In this work we adopt a mix of the above mentioned 
methods. We choose successive inclusion over successive re- 
moval of terms. The main reason for this choice is that the 
normal equations of the linear least square routine run a 
high risk of becoming degenerate when terms that respond 
similarly to the data are combined. By including terms we 
can better control the degeneracy of the normal equations. 
If degenerate normal equations were reached after the in- 
clusion of a new term, this new term was discarded since a 
possible lower order term already responded to the data in 
a similar fashion. 

Finally, we determined the local fitting functions 
through an error weighted linear least square routine (for 
individual index errors see Section I2.1.ip . Terms were suc- 
cessively inc l uded following the procedure described in 
I Gorgas et al.l (|l999[ ). by starting with the constant {j, k, I — 

in Eq. [T]) and then increasing the sum of powers j + fc + i up 
to a maximum oi j -\- k + I = ?>, including all possible cross 
terms. However, since the effective temperature is the pa- 
rameter showing the most complex behaviour we included 
polynomials of up to j = 5. If the variance was not re- 
duced at the inclusion of a new term the term was discarded. 
When a reduced variance was found the new term and all 
the previously included terms were tested by means of a 
T-test to determine if the coefficients Pi were statistically 
different from zero (by using the coefficient errors following 

1 Gorgas et al. I I 19991 and [C enarro et al.ll2002l ). Terms with co- 
efficients having a significance level a 0.1 was kept. We 
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Figure 2. The top panels a nd the lower left panel show the relationship between the chosen subregions and the analogous in the 
stellar population models of iMarastonI | |2005| ). The three panels show models for different ages (as indicated) and for each age 
the models are presented for varying metallicities as indicated by the corresponding colors. Each single model is divided into the 
main-sequence (dashed lines), sub-giant branch (solid line) and the rest of the post-main-sequence phases (dash-dotted lines). The 
lower right panel shows the relationship between the chosen subregions and the distribution of data points as a function of 6/Tgfj 
for Dwarfs (logg > 3.6, solid lined histogram) and Giants (logg > 4.0, dashed line histogram) separately. In all panels the S/Tefp 
subregions for the Dwarf subregion (solid lines) are indicated with Dl-4 and for the Giant subregion (shaded area and dashed 
lines) with Gl-4. 



then interactively studied the fitting functions and removed 
terms inducing unphysical behaviours or not affecting the 
rms scatter significantly. At the end of each run the sample 
was (j-clipped, by removing data points deviating more than 
3 a, and the fitting redone on the new sample. 

Extreme outliers that clearly deviated from the bulk of 
data points were discarded prior to running the fitting rou- 



tine. Hence to avoid stars with anomalous index strengths 
affecting the fitting functions. 

3.2 Definition of subregions in parameter space 

Thanks to the good coverage of stellar parameters the 
MILES library show a complex behaviour of the relation- 
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Table 2. Fe5335 fitting function coefficients for Lick/IDS resolution 
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Table 3. Fe5335 fitting function coefficient errors for Lick/IDS resolution 
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ship between the Lick indices and the stellar parameters. 
We have therefore divided parameter space into several sub- 
regions. 

The relationship between the Lick indices and the 
stellar parameters show a bimodality between high and 
low gravity stars (i.e. Giants and Dwarfs). The first ma- 
jor subregions that we have chosen are therefore in high 
and low values of logg space (from now on referred to as 
the Dwarf and Giant subregi o n, respectively) , in a ccor- 
dance with ' Gorgas et all l| 19991 '). ICenarro et al.l (|2002l ') and 
[Schiavon (20q3). The same log g subregion limits have been 
used for all indices. The lower limit for the Dwarf subre- 
gion was set to logg = 3.6, while the upper limit for the 
Giant subregion was set to logg = 4.0, giving an overlap 
region of A logg — 0.4. In Fig [2] the subregions are shown 
togeth er with th e analo gous in the stellar population mod- 
els of iMarastonI (|2005D . for log Q as a function of 6. The 
different evolutionary phases for the models are indicated in 
Fig [2] This shows that the choice of limits for the Dwarf and 
Giant subregions coincides very well with the division into 
the main-sequence and the post- main-sequence, as the logg 
overlap region mainly covers the sub-giant branch (SGB). 

To fully recover the detailed behaviour within the log g 



subregions we divided the full O/Tcs range into four subre- 
gions. The choice of the limits for ^/TcH subregions follow 
the behaviour of the models and the distribution of stars as 
a function O/Tss. This can be seen in Fig[2]where the limits 
of the S/Toff subregions are represented by the midpoints 
in the overlap regions, averaged over all indices. The 6/TcS 
subregions are discussed in the following bullet points, by re- 
ferring to the ^/Tcft subregions using the names (Dl-4 and 
Gl-4) in FigO first for the Giant subregion 

• Only the tip of the RGB for high metallicities fall within 
Gl (Fig[2]). The lower limit (in 8) for this subregion coincides 
with the strong drop-off in the distribution of data points 
(Fig [21 . With the weak dependency on metallicity for this 
subregion and the low number of data points we fit this 
subregion independently of metallicity. 

• G2 and G3 clearly separates out the red-giant branch 
(RGB) to be fitted mainly in G2 (Fig[2l). 

• Most indices show a distinct change in the behaviour of 
the index strengths as a function of the stellax parameters 
for hot A-type stars, around 9 = 0.5 — 0.6, see Fig. |3] and 
Fig. A1-A24. The overlap regions between G3 and G4 are 
therefore located around this range in 9. 
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Figure 3. The fitting functions for Fe5335 (Lick resolution) are shown in the upper panels for various metallicitics and over-plotted 
on data in corresponding metallicity bins. The error bars on the data are observational index errors (see Section l2.1.1l l. The colors 
correspond to [Fe/H]= -2.0 (green) , [Fe/H]= -1.35 (blue), [Fe/H]= -0.35 (red) and [Fe/H]= 0.35 (cyan) for the fitting functions 
and [Fe/H]< -1.8 (green) , -1.8 <[Fe/H]< -1.0 (blue), -1.0 <[Fe/H]< 0.2 (red) and [Fe/H]> 0.2 (cyan) for the data. The left 
and right upper panels show Giants (logg < 3.6) and Dwarfs (logg > 3.6) , respectively, for the average logg of the data in bins of 
A8 = 0.1 at steps oi 8 = 0.01. Fixed logg values are used at the ends of the 6/TcS range, with logg = 1.0, 2.0 (cold, warm end) and 
log g = 4.6, 4.0 (cold, warm end) for Giants and Dwarfs, respectively. Data points with black crosses have been sigma clipped by the 
least-square fitting routine. The lower left panel shows the residuals between the data and the fitting functions as a function of 6 
and the dashed lines represents the overall rms value for the fitting functions. The lower right panel shows the distribution of the 
residuals for three B/T^g bins, indicated by different colors where blue have 8 < 0.841 black 0.841 ^ 6 < 1.045 and red 6 ^ 1.045. 



and then for the Dwarf subregion 

• The lowest part of the main-sequence fall within Dl 
(Fig. As for the Giant subregion, the lower limit (in 9) 
for this subregion coincides with the strong drop-off in the 
distribution of data points. 

• The division of 6/TeS space into D2 and D3 were found 
to improve the fits in terms of a significantly reduced rms 
scatter. 



• As for the Giant subregion, most indices show a distinct 
change in the behaviour of the index strengths as a function 
of the stellar parameters for hot A-type stars, around 9 = 
0.5 — 0.6, see Fig. [3] and Fig. A1-A24. The overlap regions 
between D3 and D4 are therefore located around this range 
in 9. 

The number of the 9 /Tea subregions is the same for 
all indices. With the exceptions for TiOi and Ti02 that 
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show a much simpler behaviour and we have therefore used 
less e/Teff subregions (see Fig. A23-A24, Table A23-A24 and 
Table B24-B25). Since the different indices show a varying 
dependence on the stellar parameters, the limits for the sub- 
regions have been adjusted for each index individually to 
reduce the rms scatter. 

The choice of subregions in log g and O/Tcn space make 
up the base for our fitting functions. On top of these, metal- 
licity space had to be divided into two subregions for 10 
indices (CNi, CNa, Ca4227, G4300, Fe4383, Fe5015, Mgi, 
Mg2, Mgb and NaD) in order to fully reproduce the metal- 
poor end, but only in the low gravity subregion and in the 
specific temperature range around 1.0 < 9 < 1.4 (5040 < 
Toff < 3600). We have therefore independently fitted metal- 
rich and metal-poor stars, divided at [Fe/H]'--^ —1.0 for the 
affected temperatures in the low gravity subregion for the 
12 indices. 

Even though the MILES library covers an extensive 
range of stellar parameter space, the very ends are obvi- 
ously still sparsely populated. Therefore, the fitting func- 
tions are not valid beyond 6 > 1.8 (T^ff < 2800) and 6 < 0.2 
(Tcff > 25200). The dwarf main-sequence that extends to 
very low temperatures is well covered within these limits 
(Fig. [2]). Very hot young stars with temperatures greater 
than 25200 K do not have strong indices in the visual parts 
of their spectra. 



3.3 [a/Fe] trends 

Globular cluster stars are signifi cantly [q/Fej-e nhanced with 
respect to solar values (~ 0.3, ICarnevI Il996l ). The [a/Fe]- 
trend of field stars in the solar neighborhood instead show 
increasing [g/Fej-enhan cements with decre asing m etallicity 
down to [Fe/H]'^ -1.0 jEdvardsson et al. 4993: F uhrmannI 
ll998l : lMilone7 Sansom fc Sanchez-Blazauedl2009il . It is first 
at this metallicity that the field stars reach globular clus- 
ter [a/FeJ-values. Having globular cluster stars for [Fe/H]> 
— 1.0 can therefore induce [a/Fe] trends biased towards glob- 
ular cluster values in stellar libraries dominated by field 
stars. The globular cluster M71 has a metallicity of [Fe/H] = 
—0.84 and is represented by a significant number of 28 stars 
in the MILES library, which could possibly induce such a 
bias. The stars from this globular cluster were therefore dis- 
carded when computing the final fitting functions, since the 
MILES library is reasonably well populated with field stars 
around the metallicity of M71. 

The [Q/Fe]-bias of the solar neighborhood must be 
taken into account when deriving stellar population mod- 
els _based__on__empincal stellar libraries, as discussed 
in iMaraston et al. I (12003'). Model adjustments are there- 
fore needed when adopting the fitting f unctions of this 
work. S uch adjustmen ts are d escrib e d in ITripicco fc Belli 
(1 19951): iThomas et all (|2003l . 12004 ): iKorn et al l l|2005h : 
iThomas et al.l (|2005l ). 



3.4 Spectral resolution 

We have computed fitting functions for both the MILES 
and Lick/IDS resolutions (see Section [2.ip . The same final 
set of terms were used for both resolutions. Coefficients and 
coefficient errors for the fitting functions are presented in 



Appendix A for Lick resolution and Appendix B for MILES 
resolution. The sigma clipped number of data points (N) 
for the local fitting functions are also included in the coeffi- 
cient Tables, along with the rms of the residuals between the 
data and the final fitting functions, both local and overall. 
The visual behaviours, residuals and distribution of residu- 
als of the fitting functions are shown for Lick resolution in 
Appendix A. An example is presented for Fe5335 and Lick 
resolution in Table [5] and Table [3] for coefficients and coeffi- 
cient errors, respectively. The visual behaviours of the fitting 
functions for Fe5335 are shown in Fig. (3] where they are pre- 
sented for the Dwarf and Giant subregions separately and 
for varying metallicity. In Appendix A the visual behaviour 
of fitting functions for several logg values at fixed 6 axe 
also presented for indices showing strong log g dependencies 
within the logg subregions. 

3.5 Errors 

In this section we briefly discuss possible error sources af- 
fecting the final fitting functions. Such error sources include 
the index measurements of the MILES spectra, but these 
show very high quality, in terms of typical observational in- 
dex errors, as discussed in Section[2TT7l] However, the overall 
rms of the final fitting functions (see Section 13. 4p are con- 
siderably larger than the typical observational index errors 
(see Section [27lTT]). Possible error sources for this scatter are 
instead uncertainties in the stellar parameter estimates and 
intrinsic scatter in the index strengths. 

The residuals between the final fitting functions and the 
data, presented in the lower left panels of Fig. |3] and Fig. 
A1-A24 as a function of 6, show typically larger scatter for 
cooler temperatures where index values exhibit strong sen- 
sitivities to Teff. The source of this correlation is probably, 
at least partly, uncertainties in the stellar parameters, since 
these will have a larger effect when the index strengths show 
strong dependencies on the stellar parameters, i.e. 9/Tcb 
uncertainties will have less effect when the index strengths 
show weaker dependencies on O/TcS- 



4 COMPARISONS WITH THE LITERATURE 

In this section we compare the fitting functions derived in 
this work with fitting functions in the literature derived 
for stellar libraries other than MILES. We search for differ- 
ences in various parameter regimes. Comparisons are made 
wi th the classical and ex tensi vely adopted fitting fu n ctions 
of IWorthev etal] (|l994l ) and IWorthev fc Ottavianil (|l997l ) 
(from now on WFF), shifted with the offsets derived in 
Section 12.1.21 a nd with the more recent fitting functions of 
ISchiavonI (|2003) (from now on S FF) which were based on 
the JONES library (|jonej| 19991 ). 

We have performed the comparisons in different regions 
of parameter space to find the regimes where major differ- 
ences roam. The comparisons have been divided into three 
6/Tcs bins, referred to as Cold, Intermediate and Warm 
temperatures, with 9 /T^s limits presented in Table [4| Each 
of these bins have been further divided into two logy bins 
with logg = 2.0 (referred to as Giants) and logg = 4.5 (re- 
ferred to as Dwarfs) to make up a total of six bins. The aver- 
age residuals between the fitting functions were computed in 
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Table 4. Limits for the different bins of d/T^fi space used in the 
fitting function compa risons. SFF-G a nd SFF-D correspond to 
the fitting functions of ISchiavonI l|2007l ) for Giants and Dwarfs, 
respectively. 



FF COLD (Teff) INTER. (Tcff ) WARM (T^n) 

WFF 2800-4582 4582-7200 7200-12263 

SFF-G 2800-4582 4582-(~)6300 

SFF-D 3220-4582 4582-7200 7200-18000 



the metal-rich end of Cold and Warm Dwarfs that show 
stronger indices for this work beyond the 2rms level. 
SFF comparison (Fig. This work shows in general 
stronger indices in the metal-poor regime, beyond the 1 rms 
level for Intermediate temperature and Warm Dwarfs and 
beyond the 2rms level for Intermediate temperature Giants. 
In the metal-rich regime we instead find weaker for this 
work, out to the 2rms level for Intermediate temperature 
and Cold Ciants. 



each bin at [Fe/H] steps of 0.5 in the range —2 <[Fe/H]^ 0.5 
and presented in Fig.Uand Fig.[S]as a function of metallicity 
for the comparisons with WFF and SFF, respectively. 

The comparisons have only been made within the pa- 
rameter limits for which the fitting functions are applicable, 
described inLWorthey ct al. (1994) (WFF), Schiavon (2007) 
(SFF) and Section [3.21 (this work), resulting in the limits 
of the 0/Teft bins presented in Tabled Due to the limita- 
tions of the SFF we can not make comparisons for the Warm 
Giant regime, while the Intermediate Giant regime have a 
varying lower 6 limit (see .Schiavon , 2007, for individual index 
limits). 

The overall rms of the final fitting functions (see Sec- 
tion are shown in Fig. 3] -[S] as grey shaded areas (Irms 
dark grey and 2rms light grey). This gives a reference to the 
differences found between the libraries. 

Overall there is good agreement between fitting func- 
tions within the rms. We find the biggest residuals to occur 
at the ends of parameter space, i.e. at the metallicity and 
temperature ends (see Fig. U - [S]). This was expected since 
the number of data points decrease towards the ends of 
parameter space, resulting in larger uncertainties of the 
fitting functions. In the rest of this Section we discuss 
the comparisons for individual indices in terms of stellar 
parameter regions that show differences beyond the Irms 
and 2rms levels. 

WFF comparison (Fig. Warm Giants extend well 

beyond the 2rms level where this work show much weaker 
indices. We find both Warm and Cold Dwarfs to show 
stronger indices for this work, even extending beyond 
the 2rms level for the metal-poor and metal-rich ends, 
respectively. Otherwise, this work show slightly weaker 
indices extending to the Irms level. 

SFF comparison (Fig. Cold Dwarfs show weaker 

indices for this work, beyond the Irms level. Warm and 
Intermediate temperature Dwarfs show stronger indices for 
this work out to the 2rms level in the metal-poor regime. 
Intermediate temperature Giants show stronger indices 
out to the 2rms level at the ends of the metallicity scale. 
Otherwise are mainly differences within the Irms level 
found. 

H(5f 

WFF comparison (Fig. The most obvious difference is 
found for Warm Giants where this work show much weaker 
indices, extending well beyond the 2rms level. Otherwise 
are differences mainly within the Irms level, except for 



CNi 

WFF comparison (Fig. The Warm end for Giants show 
significantly stronger indices for this work, extending well 
beyond the 2rms level. Otherwise, this work show in general 
stronger indices at the metal-poor end and weaker indices 
at the metal-rich end, out to the 2rms level in both cases. 
SFF comparison (Fig.\Bj): Intermediate temperature Giants 
and Dwarfs show weaker and stronger indices for this 
work, respectively, at the metal-rich end. Otherwise are 
agreements within the Irms level mainly found. 

CN2 

WFF comparison (Fig. Similar to the previous index, 
but Warm Dwarfs show weaker indices for the entire 
metallicity scale for this work, out to the 2rms level. 
SFF comparison (Fig. Due to problems with imple- 
menting the SFFs we can not make a reliable comparison. 

Ca4227 

WFF comparison (Fig. Cold Dwarfs show stronger 
indices for this work, extending out to the 2rms level in the 
metal-poor regime, while Cold Giants instead show weaker 
indices for this work out to the 2rms level at the metal-poor 
end. Warm Giants show stronger indices for this work 
beyond the Irms level at the metal-poor end. Cold Dwarfs 
show stronger indices for this work, even extending beyond 
the 2rms level at the metal-poor end. 

SFF comparison (Fig.\Dj): The most prominent difference is 
found for Cold Giants in the metal-poor regime, extending 
well beyond the 2rms level. Cold Dwarfs show stronger 
indices for this work at the metal-rich end, beyond the Irms 
level. Otherwise are differences within the Irms level found. 

G4300 

WFF comparison (Ftg. j^)-' Warm Giants extend well 
beyond the 2rms level with stronger indices for this work. 
Metal-poor Cold Giants show stronger indices for this 
work, extending to the 2rms level. Cold metal-poor Dwarfs 
extend beyond the Irms level, showing weaker indices for 
this work. Cold and Warm metal-rich Dwarfs show stronger 
and weaker indices for this work, respectively, beyond the 
Irms level. 

SFF comparison (Fig. Differences beyond the Irms and 
2rms levels are found in several regimes, strongest at the 
ends of the metallicity scale. 

H7A 



12 J. Johansson, D. Thomas and C. Maraston 



o 

V 




o 

V 



-1 -0.5 

[Fe/H] 



0.1 

A 

cn 

* 0-05 



2 

O -0.05 




-1 -0.5 

[Fe/H] 



0-05 



5J 




-1.5 -1 -0.5 

[Fe/H] 



O -0.5 




A 



o 



o 
o 




[Fe/H] 



A 




-1.5 -1 -0.5 

[Fe/H] 



-1 -0.5 

[Fe/H] 



A 



o 



CO 
CD ■ 
CO 




-1.5 -1 -0.5 

[Fe/H] 



o 



in 

LO 

^ -0.5 

to 





o 



CD 
CO 
CD 




[Fe/H] 



[Fe/H] 



[Fe/H] 



Figure 4. C omparisons between the fitting functions of tfiis work (referred to as J09) and lWorthev et~al ] lIlQQ^ I {W94), and with 
IWorthev fc'o ttaviani (1997) (W97) for RSa, HSp, ^lA and R-yp. The panels show the difTerence J09-W94/W97 as a function of 
metaUicity for each Lick index. The comparisons are made for Giants (logg = 2.0, sohd lines) and Dwarfs (logg = 4.5, dashed 
lines). The different colors correspond to the different bins of space, with limits stated in Tablel4l where the average difference 

has been computed, blue for the Warm, black for the Intermediate and red for the Cold temperature bin. Fitting function residuals 
in terms of Irms (dark grey shaded areas) and 2rms levels (light grey shaded areas) are indicated. The errors are represented by 
the combined errors of the MILES and Lick/IDS libraries in quadrature (for more on the errors see Section EXTJ. 



WFF comparison (Fig.^: Warm Giants show significantly 
weaker indices for tliis work, well beyond the 2rms level. 
Warm Dwarfs and Cold Giants show stronger indices for 
this work out to the 2rms level. Otherwise are mainly 
differences within the Irms level found. 
SFF comparison (Fig. [S]); The most significant difference 
is found for Intermediate temperature Dwarfs, showing 
weaker indices for this work in the metal-poor regime well 
beyond the 2rms level. Otherwise are differences mainly 



found around the Irms level. 
H7F 

WFF comparison (Fig. j^).' Weaker indices are found for 
this work for Warm Giants and Cold Dwarfs beyond the 
2rms and Irms level, respectively. Otherwise are mainly 
differences within the Irms level found. 
SFF comparison (Fig. Due to problems with imple- 
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Figure 5. Same as Fig. [4] but for the comparison between this work (J09) and lSchiavoiil ll2007t) (S07). The errors are represented 
by index errors of the MILES library (see Section 12.1.111 . 



meriting the SFFs we can not make a reliable comparison. 
Fe4383 

WFF comparison (Fig. Warm Giants and Warm 

Dwarfs show stronger indices out to the 2rms level at the 
metal-poor end. Cold Giants show weaker indices, out 
to the 2rms level at the metal-poor and metal-rich ends. 
Otherwise are mainly diflterences within the Irms level 
found. 

SFF comparison (Fig. Cold Dwarfs show significantly 
stronger indices for this work, well beyond the 2rms level. 
Warm Dwarfs instead show weaker indices, out to the 2rms 
level. Otherwise are mainly differences within the Irms 
found. 

Ca4455 

WFF comparison (Fig. Cold Giants show stronger 
indices for this work, extending well beyond the 2rms level 
in the metal-rich regime. Cold Dwarfs show weaker indices 
for this work, extending well beyond the 2rms level. The 
Warm regime show stronger indices for this work, extending 
beyond the Irms level in the metal-poor regime. 



Fe4531 

WFF comparison (Fig. Cold Dwarfs show stronger 
indices for this work at the Irms level for the metal-poor 
end and increasing well beyond the 2rms level at the 
metal-rich end. Intermediate temperature Giants show 
stronger indices for this work, out to the 2rms level at the 
metal-rich end. 

C24668 

WFF comparison (Fig.^: The metal-rich end show weaker 
indices for this work, mainly down to the Irms level. The 
metal-poor end show weaker and stronger indices for this 
work extending beyond the Irms level for the Warm bins 
and Cold Dwarfs, respectively. 

SFF comparison (Fig.\^: Cold Giants show weaker indices 
for this work beyond the Irms level in the metal-rich 
regime. Intermediate temperature Dwarfs show stronger 
indices for this work at the metal-rich end out to the 2rms 
level. 

WFF comparison (Fig. We find this work to show 



New flux- calibrated Lick fitting functions 15 




[Fe/H] [Fe/H] [Fe/H] 




Figure 5 - continued 



weaker indices for Warm Giants well beyond the 2rms 
level. Cold Dwarfs show stronger indices for this work 
beyond the 2rms level. Cold Ctants show weaker indices for 
this work beyond the Irms level. 

SFF comparison (Ftg. Cold Dwarfs show the biggest 
differences right beyond the 1 rms level. Otherwise are 
differences within the Irms level are found. 

Fe5015 

WFF comparison (Fig.^: Cold Giants show weaker indices 
for this work, beyond the 2rms level in the metal-poor 
regime. Intermediate temperature Giants show stronger 
indices for this work beyond the 2rms level in the metal-rich 
regime. 

SFF comparison (Fig. Cold Giants show significantly 
weaker indices for this work, well beyond the 2rms level. 
This work shows weaker indices for Cold Dwarfs beyond 
the Irms level in the metal-poor regime. Otherwise are 
differences within the Irms mainly found. 

Mgi 

WFF comparison (Fig.^: The Cold Dwarfs show stronger 
indices for this work, out to the 2rms level in the metal-poor 
regime. The Cold Giants show stronger indices for this 
work at intermediate metallicities. No differences found 
beyond the 2rms level. 



SFF comparison (Fig. 0).' The Cold end show weaker 
indices for this work in the metal-poor regime, beyond the 
2rms level. Cold Dwarfs instead show stronger indices for 
this work beyond the 2rms level at the metal-rich end. 

Mgb 

WFF comparison (Fig. Cold and Intermediate temper- 
ature Giants show weaker indices for this work, extending 
beyond the 2rms level. Intermediate temperature and Cold 
Dwarfs show weaker and stronger indices, respectively, for 
this work in the Metal-poor regime, beyond the Irms level. 
SFF comparison (Fig. [3P; Due to problems implementing 
the SFFs we can not make a reliable comparison. 

Fe5270 

WFF comparison (Fig. The Warm end show stronger 
indices beyond the 2rms level in the metal-poor regime. 
Cold Giants show weaker indices for this work beyond the 
2rms level. The metal-rich end show stronger indices for this 
work beyond the Irms level for intermediate temperature 
cold Dwarfs. 

SFF comparison (Fig. [3P; We find weaker indices for this 
work beyond the Irms level for Cold Giants. Stronger 
indices for this work beyond the 2rms level are found 
for Cold Dwarfs in the intermediate metallicity regime. 
Otherwise are differences well within the Irms level found. 



Mg2 

WFF comparison (Fig. Cold Dwarfs show stronger 
indices for this work beyond the Irms level in the metal- 
poor regime. Cold Giants show stronger indices for this 
work beyond the Irms level for intermediate metallicities. 
Otherwise are mainly differences within the Irms level 
found. 



Fe5335 

WFF comparison (Fig.^: Cold Giants show weaker indices 
for this work well beyond the 2rms level. The Warm end 
show stronger indices out to the Irms level at intermediate 
metallicities. Cold Dwarfs show weaker and stronger indices 
out to the Irms level a the metal-poor and metal-rich ends, 
respectively. 
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SFF comparison (Fig. Cold Dwarfs show stronger 

indices beyond the 2rms level at intermediate metallicities 
and weaker indices beyond the 2rms level at the metal-rich 
end. Cold Gtants show weaker indices for this work beyond 
the Irms level. Intermediate temperature Dwarfs instead 
show stronger indices for this work, beyond the Irms level 
in the metal-rich regime. 

Fe5406 

WFF comparison (Fig. The Cold end show stronger 
indices for this work beyond the Irms level in the metal-rich 
regime. Intermediate temperature Giants show weaker 
indices for this work beyond the Irms level. No differences 
found beyond the 2rms level. 

Fe5T09 

WFF comparison (Fig. Cold Dwarfs show stronger 
indices for this work, extending beyond the 2rms level at 
the metal-rich end. Otherwise, no significant differences 
beyond the Irms level. 

Fe5782 

WFF comparison (Fig. Warm Giants show stronger 
indices for this work in the metal-poor regime, extending 
beyond the 2rms level at the metal-poor end. Cold Giants 
show stronger indices for this work regime, beyond the Irms 
level. Warm Dwarfs show stronger indices for this work 
beyond the Irms level at the metal-poor end. 

NaD 

WFF comparison (Fig. The Cold end show stronger 
indices for this work extending well beyond the 2rms 
level, especially in the metal-poor regime. Intermediate 
temperature Giants show stronger indices for this work, 
extending beyond the 2rms level in the metal-rich regime. 

TiOi 

WFF comparison (Fig. Cold Dwarfs show weaker 

indices for this work, extending well beyond the Irms level 
at the metal-rich end. Otherwise, no differences found 
beyond the Irms level. 

TiOa 

WFF comparison (Fig. We find significantly weaker 
indices for this work for Cold Dwarfs, extending very far 
beyond the 2rms level. Otherwise, no significant differences 
found beyond the Irms level. 



the MILES library and for the resolution of the Lick/IDS 
library. The MILES library consists of 985 stars selected 
to produce a sample with extensive stellar parameter cov- 
erage. The MILES library was also chosen because it has 
been carefully flux-calibrated, making standard star derived 
offsets unnecessary. This becomes important when compar- 
ing stellar population models to high redshift data where no 
resolved individual stars are available. 

We find the index measurements of the MILES spec- 
tra to have very high quality in terms of observational in- 
dex errors. These errors are also found to be significantly 
smaller than for the Lick/IDS library. This was expected 
since the MILES library was observed nearly thirty years 
after the Lick/IDS library. Given the high quality of the 
index measurements, index errors should not be the major 
error sources for the final fitting functions. We instead find 
indications that the stellar parameter estimates are signifi- 
cant error sources. 

Lick Index offsets between the MILES library and the 
classic Lick/IDS library are derived in order to be able to 
compare stellar population models based on this work with 
models in the literature. We find these offsets to be depen- 
dent on index strength and have therefore derived least- 
square fits for the residual between the two libraries. Off- 
set to the Lick/IDS library are also derived for the flux- 
calibrated ELODIE and STELIB libraries. We flnd clear off- 
set deviations between the libraries. The largest deviations 
are found for the STELIB library compared to the other 
two libraries, which is also the library having least stars in 
common with the Lick/IDS library. The deviations in offsets 
found between the three libraries undermine the derivation 
of universal offsets between the Lick/IDS and these flux- 
calibrated systems. 

We compare the fitting functions of this work to fit- 
ti ng functions in the lite r ature, namely the fitting func tions 



t mg tunctions m tne ute r ature, namely tne nttmg tunc tions 
of IWorthev et al.1 l| 19941 ). IWorthev fc Ottavianil l|l997l ) and 
ISchiavonI l|2007l ). Generally we find good agreement within 
the rms of the residuals between the data and the fitting 
functions of this work. The differences found in the compar- 
isons vary significantly from index to index and especially 
from one stellar parameter region to another for individual 
indices. However, the major differences are found in the out- 
skirts of stellar parameter space, i.e. at the temperature and 
metallicity ends. This is probably due to a low number of 
data points in these regimes for the stellar libraries, inducing 
uncertainties which result in the major differences found. 

In a forthcoming paper (Thomas et al. in prep.) the 
fitting functions of this work will be implement ed in stel- 
lar population models following the techniques of iMarastonI 
(2005) and Thomas et al.. (,2003i l. 

A user friendly Fortran 90 code is available online at 
www.icg.port.ac.uk/~johanssj to easy the implementation 
of our fitting functions in population synthesis codes. 



5 SUMMARY 

We have derived new empirical fitting functions for the rela- 
tionship between Lick absorption indices and stellar atmo- 
spheric parameters (Teg, [Fe/H] and logg) described by the 
MILES library of stellar spectra, both for the resolution of 
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APPENDIX A: FITTING FUNCTIONS FOR LICK/IDS RESOLUTION 

Short version, full appendix can be found at www.icg.port.ac.uk/~johanssj 




Figure Al. H15A, same as Fig. 3. Dashed-lines for additional log g- value (log g = 3.3) cover the data points in strong log g-dependent 
regions. 
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APPENDIX B: FITTING FUNCTIONS FOR MILES RESOLUTION 

Short version, full appendix can be found at www.icg.port.ac.uk/~johanssj 
Table Bl. HSa fitting function coefficients for MILES resolution 

overall rms=0.8009 

logg ^ 4.0 and 9 limits: logf; ^ 3.6 and 9 limits: 



Term 


0.62 


0.53 - 1.1 


0.90 - 1.4 


^ 1.3 


^ 0.57 


0.50 - 0.90 


0.80 - 1.2 


^ 1.08 


Const. 


-8.002 


-136.8 


36.39 


-127.4 


-3.369 


-134.7 


-49.75 


-41.96 


6 


-12. 19 


487.8 


X 


23.06 


16.05 


596.3 


297.8 


42.42 


[Fe/H] 


X 


-2.367 


41.22 


X 


X 


-13.31 


31.61 


X 


log 9 


-4.330 


21.43 


2.662 


X 


X 


X 


-0.5889 


X 


e= 


X 


-550.1 


-88.27 


103.1 


13.26 


-922.4 


-423.0 


-9.752 


[Fe/H] 2 


X 


X 


0.5080 


X 


X 


-3.267 


X 


X 




5.072 


X 


X 


X 


0.1957 


X 


X 


X 


9 [Fe/H] 


X 


X 


-78.94 


X 


X 


31.04 


-68.15 


X 


log 9 


X 


-44.24 


-6.149 


X 


0.2840 


26.90 


X 


X 


[Fc/H]log9 


X 


0.2318 


X 


X 


X 


X 


X 


X 


93 


273.0 


195.6 


48.49 


49.23 


X 


448.91 


175.3 


X 


[Fc/H]3 


X 


0.5286 


0.3213 


X 


X 


-0.1733 


-0.1461 


X 




-0.8813 


X 


X 


X 


X 


X 


X 


X 


e^[Fe/H] 


X 


X 


34.99 


X 


X 


-20.90 


33.01 


X 


9^ log g 


-16.26 


22.55 


3.167 


X 


X 


-17.20 


-1.042 


X 


9 [Fe/H] 2 


X 


1.297 


X 


X 


X 


2.919 


-0.7101 


X 


[Fe/H] 2 log B 


X 


X 


X 


X 


X 


0.02043 


X 


X 


91ogg^ 


2.150 


X 


X 


X 


X 


-1.628 


X 


X 


9* 


-275.2 


X 


X 


-108.1 


X 


X 


X 


X 


9= 


X 


X 


X 


31.52 


X 


X 


X 


X 


rms 


1.033 


1.076 


0.8149 


0.6510 


0.8884 


0.5558 


0.5954 


1.477 


N 


96 


346 


359 


41 


49 


278 


276 


33 



